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SUMMARY 


A  model  of  a  simple  paral 1  el -shaft ,  spur-gear  transmission  is  presented. 


.1  ...vs  N.  •  '-'I  ■»  '  '  »»  I  I  V  ■  j  I  ’  VA  1  V.  ,  ^  V-  VA  I  1.1  nil  JIM  I  J  J  iwil  1  J  piV.Jtlll.tXJ. 

The  model  is  developed  to  simulate  dynamic  loads  in  power  transmissions.  Fac¬ 
tors  affecting  these  loads  are  identified.  Included  are  shaft  stiffness  and 
inertia,  load  and  power  source  inertia,  tooth  geometry,  tooth  stiffness,  local 
compliance  due  to  contact  stress,  load  sharing,  and  friction.  Governing  dif¬ 
ferential  equations  are  developed  and  a  solution  procedure  is  outlined.  A 
parameter  study  of  the  solutions  is  presented  in  NASA  TM-100181  (AVSCOft 
TM-87-C-3) .  '\ 


INTRODUCTION 


Recently,  there  has  been  increased  interest  in  the  dynamic  effects  in 
gear  systems.  This  interest  is  stimulated  by  demands  for  stronger,  higher 
speed,  improved-performance,  and  longer  lived  systems.  This  in  turn  has  stimu¬ 
lated  numerous  research  efforts  directed  toward  understanding  gear  dynamic  phe¬ 
nomena.  However,  many  aspects  of  gear  dynamics  are  still  not  satisfactorily 
understood. 


For  example,  in  industrial  settings,  a  high  performance  gear  system  is 
often  obtained  by  overdesigning,  and  by  sacrificing  costs,  materials,  and  com¬ 
pactness.  In  aerospace  and  military  application  where  weight  is  a  premium, 
gear  systems  are  often  designed  under  conditions  very  close  to  the  failure  lim¬ 
its,  thereby  introducing  uncertainties  in  performance  and  life  prediction. 

They  are  often  prematurely  replaced  to  prevent  in-service  failure.  Moreover, 
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gear  systems  are  often  designed  by  using  static  analyses.  However,  when  gear 
systems  operate  at  high  speed,  there  are  several  factors  which  affect  their 
performance.  These  include  shaft  torsional  stiffness,  gear  tooth  loading  and 
deformation,  gear  tooth  spacing  and  profile  errors,  rotating  speeds,  mounting 
alignment,  dynamic  balance  of  rotating  elements,  gear  and  shaft  masses  and 
inertia,  and  the  masses  and  inertias  of  the  driving  (power)  and  driven  (load) 
elements . 

There  is  no  agreement  among  researchers  on  the  best  methods  for  evalu¬ 
ating  dynamic  load  effects.  Hence,  gear  designers  are  often  confronted  with 
conflicting  theories.  They  generally  have  to  rely  on  past  experience,  service 
safety  factors,  and  experimental  data  with  a  limited  range  of  applicability. 

The  objective  of  this  report  and  NASA  TM-100181  (AVSCOM  TM-87-C-3)  is  to 
provide  more  insight  into  the  factors  affecting  dynamic  loads. 

Research  efforts  on  gear  system  dynamics  have  been  conducted  for  many 
years.  In  1892  Lewis  (ref.  1)  recognized  that  the  instantaneous  tooth  load 
was  affected  by  the  velocity  of  the  system.  In  1925,  Earl  Buckingham  (ref.  2) 
headed  an  experimental  research  effort,  endorsed  by  ASME,  to  measure  dynamic 
effects.  A  report  published  in  1931  represented  the  first  authoritative  docu¬ 
ment  on  gear  dynamics.  It  presented  a  procedure  for  determining  the  so-called 
dynamic  load  increment  due  to  mesh  dynamics  and  gear  tooth  errors. 

In  1959,  Attia  (ref.  3)  performed  an  experiment  to  determine  actual 
instantaneous  loading.  He  found  that  Buckingham's  results  were  conservative. 

In  1958,  Niemann  and  Rett i g  (ref.  4)  found  that  larger  masses  caused 
higher  dynamic  loads,  but  as  the  average  load  became  larger  the  effect  of 
larger  masses  became  less  important.  They  also  found  that  very  heavily  loaded 
gear  systems  showed  no  appreciable  dynamic  load  increment,  whereas  in  lightly 
and  moderately  loaded  gear  systems  there  were  considerable  dynamic  load  incre¬ 
ments.  In  1958,  Harris  (ref.  5)  suggested  that  for  gear  systems  isolated  from 
external  stimuli,  there  are  three  internal  sources  of  dynamic  loads  as  follows 

(1)  error  in  the  velocity  ratio  measured  under  the  working  load 

(2)  parametric  excitation  due  to  stiffness  variation  of  the  gear  teeth 

(3)  nonlinearity  of  tooth  stiffness  when  contact  is  lost 

In  1970,  Houser  and  Seireg  (ref.  6)  developed  a  generalized  dynamic  factor  for 
mula  for  spur  and  helical  gears  operating  away  from  system  resonances.  The 
formula  took  into  consideration  the  gear  geometry  and  manufacturing  parameters 
as  well  as  the  dynamic  characteristics  of  the  system. 

In  1972,  Ichimaru  and  Hirano  (ref.  7)  analyzed  heavy-loaded  spur-gear  sys 
terns  with  manufacturing  errors  under  different  operating  conditions.  They 
found  that  the  change  in  tooth  profile  showed  a  characteristic  trend  to 
decrease  dynamic  load.  In  1978,  Cornell  and  Westervelt  (ref.  8)  presented  a 
closed  form  solution  for  a  dynamic  model  of  spur-gear  system  and  showed  that 
tooth  profile  modification,  system  inertia  and  damping,  and  system  critical 
speeds,  can  have  significant  effects  upon  the  dynamic  loads.  In  1981,  Kasuba 
and  Evans  (ref.  9)  presented  a  large  scale  digitized  extended  gear  modeling 
procedure  to  analyze  spur-gear  systems  for  both  static  and  dynamic  conditions. 
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Their  results  indicated  that  gear  mesh  stiffness  is  probably  the  key  element 
in  the  analysis  of  gear  train  dynamics.  They  showed  that  the  gears  and  the 
adjacent  drive  and  load  systems  can  be  designed  for  optimum  performance  in 
terms  of  minimum  allowable  dynamic  loads,  for  a  wide  range  of  operating  speeds. 

In  1981,  Wang  and  Cheng  (ref.  10)  developed  another  dynamic  load  response 
algorithm.  They  reported  that  the  dynamic  load  is  highly  dependent  on  the 
operating  speed.  This  model  is  later  modified  by  Lewicki  (ref.  11)  to  account 
for  the  nonlinear  Hertzian  deformation  of  meshing  gear  teeth.  The  gear  dynamic 
load  found  from  the  revised  model  showed  little  difference  from  the  original 
model  since  the  Hertzian  deflection  was  relatively  small  in  comparison  with 
the  total  gear  tooth  deflection.  Nagaya  and  Uematsu  (ref.  12)  stated  that 
because  the  contact  point  moves  along  the  involute  profile,  the  dynamic 
response  should  be  considered  as  a  function  of  both  the  position  and  speed  of 
the  moving  load.  In  1982,  Terauchi ,  et  a  1 . ,  (ref.  13)  studied  the  effect  of 
tooth  profile  modifications  on  the  dynamic  load  of  spur-gear  systems.  Accord¬ 
ing  to  their  results,  the  dynamic  load  decreased  with  proper  profile 
modifications . 

In  this  first  part  of  the  paper,  we  present  a  model  of  a  parallel  shaft 
transmission.  We  consider  the  effects  of  shaft  stiffness  and  inertia,  load 
and  power  source  inertias,  tooth  stiffness,  local  compliance  due  to  contact 
stresses,  load  sharing,  and  friction.  A  parameter  study  is  provided  in  the 
second  part  of  the  paper. 


NOMENCLATURE 

cross  section  area  of  i th  element  of  gear  teeth,  mm2  (in.2) 
damping  coefficient,  gear  tooth  mesh,  N-sec  (lb-sec) 
damping  coefficient  of  shaft,  N-m-sec  (in. -lb-sec) 
effective  modulus  of  elasticity,  N/m2  (lb/in.2) 
tooth  face  width,  mm  (in.) 
shear  modulus,  N/m2  (lb/in.2) 

second  moment  of  inertia  of  i th  element  of  gear  teeth,  mm2  (in.2) 

polar  moment  of  inertia  of  load,  m2-kg 

polar  moment  of  inertia  of  motor,  m2-kg 

polar  moment  of  inertia  of  gear  1,  m2-Kg 

polar  moment  of  lertia  of  gear  2,  m2-Kg 

stiffness  of  gear  tooth,  N-m/rad 

stiffness  of  shaft,  N-m/rad 
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distance  between  elements  i  and  j ,  mm 

gear  tooth  deformation  due  to  beam  deflection,  mm 

gear  tooth  deformation  due  to  contact  deformation,  mm 

gear  tooth  deformation  due  to  foundation  flexibility,  mm 

total  gear  tooth  deformation,  mm  (in.) 

base  radii  of  gears,  mm  (in.) 

pi tch  radi i  of  gears ,  mm 

thickness  of  element  i,  mm 

torque  on  load,  N-m 

torque  on  motor,  N-m 

torque  on  gear  1 ,  N-m 

torque  on  gear  2,  N-m 

sliding  velocity  during  tooth  mesh,  mm/sec 

rolling  velocity,  mm/sec 

applied  load,  N/m  (Ib/in.) 

x-coordinate  of  element  i,  mm  (in.) 

y-coordinate  of  element  i,  mm  (in.) 

load  angle,  rad 

backlash ,  mm  (in.) 

angular  displacement,  rad 

angular  velocity,  rad/sec 

angular  acceleration,  rad/sec- 

lubricant  viscosity,  N-sec/m2  ( 1 b-sec/ i n . 2 ) 

Poisson's  ratio 
damping  ratio 


MODELING 


Figure  1  depicts  a  model  of  the  transmission.  It  consists  of  a  motor  or 
power  source  connected  by  a  flexible  shaft  to  the  gear  system.  The  gear  sys¬ 
tem  consists  of  a  pair  of  involute  spur  gears.  They  are  connected  to  the  load 


by  a  second  flexible  shaft  as  shown.  Symbolically,  the  model  may  be  repre 
sented  by  a  collection  of  masses,  springs,  and  dampers  as  in  figure  2. 


Let  ©m,  0],  02,  and  ©[_  represent  the  rotations  of  the  motor,  the 
gears,  and  the  load.  Then  by  using  standard  procedures  of  analysis,  the  gov¬ 


erning  differential  equations  for  the  rotations  may  be  written  as 

J  +  C  s  ]  <  ©m  -  0  ]  )  •'T  s  1  f  ©M  _  ©  1  ^  =  T  M  <1) 

J 1  ©  1  +  Csi(©i  -  0^)  ♦  Ks](©i  -  0^)  +  Cg  <  t )  [  R  t)  i  ©  i  -  Rb2®2^ 

+  Kg(t)[Rpi(Rp]0]  -  Rb2©2^  =  T  f  i  ( t )  (2) 

J 2©2  +  C s 2 ^ _  ©1 )  +  K s 2 ^ ®2  _  ©1 )  +  Cg(t)[Rp2©2  _  R b 1  ©  1 ^ 

+  Kg ( t )  [ Rp2 ( Rb2®2  “  Rb  1  el >  3  =  Tf2<t)  (3) 

J[_0[_  +  CS2(9l  _  ©2 ^  +  Ks2^©L  -  ©2 )  =  ~T[_  (4) 


where  J^,  J i  ,  J2 ,  and  J|_  represent  the  mass  moments  of  inertia  of  the 
motor,  the  gears,  and  the  load;  Cs 1 ,  CS2,  and  Cg(t)  are  damping  coefficients 
of  the  shafts  and  the  gears;  Ks],  KS2,  and  Kg(t;  are  stiffnesses  of  the 
shafts  and  the  gears;  ,  T|_,  Tf](t),  and  Tf  2  <  t )  are  motor  and  load  torques 
and  frictional  torques  on  the  gears;  Rpi  and  Rb2  are  base  circle  radii  of 
the  gears;  t  is  time;  and  the  dots  over  ©  indicate  time  di f ferent i at  ion . 

In  developing  equations  (1)  to  (4)  several  simplifying  assumptions  are 
employed : 

(1)  The  dynamic  process  is  studied  in  the  rotating  plane  of  the  gears. 
Out-of-plane  twisting  and  misalignment  are  neglected. 

(2)  Damping  due  to  lubrication  of  the  gears  and  shafts  is  expressed  in 
terms  of  constant  damping  factors. 

(3)  The  differential  equations  of  motion  are  developed  by  using  the  theo¬ 
retical  line  of  action. 

(4)  Low  contact  ratio  gears  are  used  in  the  analysis.  Specifically,  the 
contact  ratio  is  taken  between  1  and  2. 


ANALYSIS 

A  major  task  in  the  analysis  is  to  determine  the  values  of  the  stiffness, 
damping,  and  friction  coefficients  appearing  in  equations  (1)  to  (4).  Another 
task  is  to  determine  the  ratio  of  load  sharing  between  the  teeth  during  a  mesh 
cycle.  These  factors  depend  on  the  roll  angle  of  the  gears.  Thus,  equations 
(1)  to  (4)  are  made  nonlinear  bv  these  terms. 


Stiffness 


Gear  stiffness.  -  Consider  first  the  stiffness  coefficients  kg,  ks], 
and  kS2-  Let  the  tooth  surface  have  the  form  of  an  involute  curve.  Let  Wj 
be  the  transmitted  load  at  a  typical  point  j  of  the  tooth  profile.  Let  qj 
be  the  deformation  of  the  tooth  at  point  j  in  the  direction  of  Wj .  Then  J 
the  gear  stiffness  kg j  for  the  gear  teeth  in  contact  at  j  is 


(5 


In  general  qj  will  depend  on  the  following:  (1)  the  bending  of  the 
tooth,  the  shear  deformation  of  the  tooth,  and  the  axial  compression  of  the 
tooth;  (2)  the  deflection  due  to  the  flexibility  of  the  tooth  foundation  and 
fillet;  and  (3)  the  local  compliance  due  to  the  contact  stresses. 


To  determine  qj  let  the  tooth  be  divided  into  elements  as  shown  in 
figure  3.  Let  i  be  a  typical  element  with  thickness  Tj,  cross  section 
area  A, ,  and  second  moment  of  inertia  Ij.  Let  Lm  be  the  distance  between 
element  i  and  point  j  along  the  x-axis.  Let  (3j  be  the  angle  between 
Wj  and  the  y-axis.  (See  fig.  3.) 

Consider  the  tooth  to  be  a  nonuniform  cantilever  beam.  Let  q^j  be  the 
contribution  to  qj  by  the  bending,  shear,  and  axial  deformation  of  the 
tooth.  Then  q^j  may  be  represented  as  the  sum  of  the  deformation  in  the  ele 
ments  i  beneath  point  j.  That  is, 


n 


i  =  1 


(6 


where  n  is  the  number  of  elements  beneath  j  and  q^ i j  is  the  deformation 
of  element  i  due  to  the  load  Wj . 


With  standard  analysis  qpn  is  found  to  have  the  value  (refs.  8,  14, 
and  15): 


q,  .  .  =  (W./E  ) 
bi  j  J  e 


cos2f3. 


T.3  +  3T.2L.  .  +  3T.L2.  .  j/31 . 
i  1  1 J  i  i  J  //  i 


cos  13 j  sin  (3 ^ 


T . “Y  .  +  2T.Y.L. .  2 1 . 

1  J  1  J  1  J  //  1 


+  cos  (3^ 


12(1+  v)T. / 5A . 

+  s i  n2|3 .  (T.  /A.  )  { 
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where  Y j  is  the  half-tooth  thickness  at  element  i  (see  fig.  3),  v  is 
Poisson's  ratio,  and  Ee  is  the  "effective  elastic  modulus"  depending  upon 
whether  the  tooth  is  wide  (plane  strain)  or  narrow  (plane  stress).  Specifi¬ 
cally,  for  a  "wide"  tooth,  where  the  ratio  of  the  width  to  thickness  at  the 
pitch  point  exceeds  5  (ref.  14),  Ee  is 


where  E  is  Young's  modulus  of  elasticity.  For  a  "narrow"  tooth  (width- to- 
thickness  ratio  less  than  5),  Ee  is 

Ee  =  E  (9: 

Expressions  similar  to  equations  (7)  hold  for  qfj  the  contribution  to 
qj  for  the  deformation  due  to  the  flexibility  of  the  tooth  fillet  and  founda¬ 
tion  (ref.  15). 

Let  qcj  be  the  contribution  to  qj  from  the  local  compliance  due  to 
contact  stresses.  With  the  procedures  or  Lundberg  and  Palmgren  (ref.  16)  qCj 
may  be  expressed  as 

acj  - 

J  rO.9rO.8n  0.1 
E  F  W . 

J 

where  F  is  the  width  of  the  tooth. 

Hence,  by  superposition,  the  deformation  at  j  in  the  direction  of  Ws 


q.  .  qbj  .  .  qcJ  (11) 

The  above  expressions  were  used  to  calculate  the  deformations  for  two  dif¬ 
ferent  gear  pairs.  The  results  are  shown  graphically  in  figures  4(a)  and  (b). 

Shaft  stiffness.  -  The  shaft  stiffness  Ks  is  given  by  the  standard 
expression 

Ks  ■  r  <’21 

where  G  is  the  shear  modulus,  Q.  is  the  shaft  length,  and  J  is  the  polar 
moment  of  area  given  by 


where  D  is  the  shaft  diameter 


Damp i ng 

Shaft  damping.  -  Next  consider  the  damping  coefficients  CS],  CS2,  and 
Cg.  Damping  in  the  shafts  is  due  to  the  shaft  material.  In  equations  (1)  to 
(4)  the  coefficients  Csi  and  CS2  are  taken  to  have  the  form 


.  j*.  j  .  .  »*  jl. 
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csi  =  2^s  ul: 

JM  +  J1 


Cs2  =  2^sLl  (15: 

(JL  +  J2) 

where  Es  represents  the  damping  ratio.  Experiments  have  shown  that  has 

values  between  0.005  and  0.075  (ref.  17). 

Mesh  damping.  -  Similarly,  the  effect  of  damping  of  the  gear  mesh  is 
taken  as 


K  R2  R2  J,J, 
f'  9  c  9  u  I  D  2  I  2 

Ln  =  ^  2  2 

9  Rblh  *  Rb2J2 


where,  as  before,  E  is  the  damping  ratio.  Measurements  have  shown  £  to 
have  values  between  0.03  and  0.17  (refs.  9  and  10). 

Friction.  -  Equations  (1)  to  (4)  contain  terms  Tfi  and  Tf2  which  rep¬ 
resent  the  frictional  moments  of  the  driving  and  driven  gears.  These  moments 
occur  because  of  the  relative  sliding  of  the  gear  teeth.  Buckingham  (ref.  18) 
has  recorded  a  semi empi r i cal  formula  for  the  friction  coefficient  f  of 
boundary  lubrication  as 


f  =  0.05e~°'  125Vs&  +  0.002 


where  Vsa  is  the  sliding  speed  measured  in  in. /sec.  An  analogous  expression 
for  elastohydrodynamic  lubrication  has  been  developed  by  Beneaict  and  Kelly 
(ref.  19)  and  by  Anderson  and  Loewenthal  (ref.  20)  as  follows: 


f  =  0.0127  log  45.94  W/Fn  V  nV" 

U  S  xt  K 


where 


W  the  applied  load,  N/m  (lb/in.) 

F  face  width,  mm  (in.) 

Vr  rolling  velocity,  mm/sec  (in. /sec) 

Po  lubricant  viscosity,  N-sec/m2  ( 1 b-sec/ i n . 2 ) 


Figures  5(a)  and  (b)  show  graphs  of  the  friction  coefficient  as  given  by 
equations  (17)  and  (18)  as  a  function  of  the  roll  angle.  Figures  6(a)  and  (b) 
show  the  resulting  effect  upon  the  friction  torque. 


Mesh  Ana  lysis 

Figure  7  illustrates  the  motion  of  a  pair  of  meshing  teeth.  The  initial 
contact  occurs  at  A,  where  the  addendum  circle  of  the  driven  gear  intersects 
the  line  of  action.  As  the  gears  rotate  the  point  of  contact  will  move  along 
the  line  of  action  APD.  When  the  tooth  pair  reaches  B,  the  recessing  tooth 
pair  disengages  at  D  leaving  only  one  zone.  When  the  tooth  pair  reaches 
point  C,  the  next  tooth  pair  begins  engagement  at  A  and  starts  another 
cycle . 

In  the  analysis,  the  position  of  the  contact  point  of  the  gear  teeth 
along  the  line  of  action  is  expressed  in  terms  of  roll  angles  of  the  driving 
gear  tooth. 

Figure  8  shows  typical  stiffness  and  load  sharing  characteristics  through 
a  mesh  cycle.  Let  a  series  of  mating  tooth  pairs  be  denoted  as  a,  b,  c,  d 
and  let  points  A,  B,  P,  C,  D  be  the  same  as  those  in  figure  7.  Then,  AB 
and  CD  represent  the  double  contact  regions,  BC  represents  the  single  con¬ 
tact  region,  and  as  before  P  is  the  pitch  point. 

The  stiffness  values  at  double  contact  regions  are  cleariy  much  higher 
than  those  at  single  contact  regions.  When  gears  rotate  at  appreciable  speed, 
this  time-varying  stiffness  as  shown  in  figure  8  is  the  major  excitation 
source  for  the  dynamic  response  of  the  system. 


DISCUSSION 

The  objective  of  this  analysis  is  to  establish  the  governing  differential 
equations  and  to  present  a  procedure  for  solution.  As  noted,  the  equations 
themselves  are  nonlinear.  However,  they  may  be  efficiently  solved  by  using  a 
1 i near i zed- i terat i ve  procedure  as  follows: 

The  linearized  equations  may  be  obtained  by  dividing  the  mesh  cycle  into 
n  equal  intervals.  Let  a  constant  input  torque  T^  be  assumed.  Let  the  out¬ 
put  torque  T|_  be  fluctuating  because  of  damping  in  the  gear  mesh,  because  of 
friction,  and  because  of  time-varying  mesh  stiffness. 

Let  initial  values  of  the  angular  displacements  be  obtained  by  preloading 
the  input  shaft  with  the  nominal  torque  carried  by  the  system.  Initial  values 
of  the  angular  speeds  may  be  taken  from  the  nominal  operating  speed  of  the 
system. 

The  iterative  process  is  then  as  follows:  the  calculated  values  of  the 
angular  displacements  and  angular  speeds  after  one  period  are  compared  with 
the  assumed  initial  values.  Unless  the  differences  between  them  are  suffi¬ 
ciently  small,  the  procedure  is  repeated  by  using  the  average  of  the  initial 
and  calculated  values  as  new  initial  values. 


Finally,  observe  that  the  term  ( R^, i ©i  -  Rb2e2}  'n  the  equation  of  motion  ' 
represents  the  relative  dynamic  displacement  of  the  gears.  Let  S  represent  \ 
the  backlash.  Let  gear  1  be  the  driving  gear.  The  following  conditions  can  : 
occur  < 


Case  1 

The  normal  operating  case  is 

Rblel  -  Rb2©2  >  0  (19) 

The  dynamic  mesh  force  F  is  then 

F  =  Kg( t) (Rbl©l  -  Rb2e2)  +  Cg(t)<Rbi©i  -  Rb2®2^  (20) 

Case  2 

Rbl®l  ~  Rb2®2  =  0  and  |  Rb l © l  -  Rb2®2 1  =  &  (21) 

In  this  case,  the  gears  will  separate  and  the  contact  between  the  gears  will 
be  lost.  Hence, 

F  =  0  (22) 

Case  3 

Rbl®l  -  Rb2®2  <  0  and  I R b 1 ® 1  ~  Rb2®2l  >  &  (23) 

In  this  case,  gear  2  will  collide  with  gear  1  on  the  backside.  Then, 

F  =  Kg( t ) [ ( Rb202  -  Rb 1 e 1 }  -  +  Cg(t )<Rbl©l  -  Rb2^2)  (24) 
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CONCLUSION 

A  low  contact  ratio  spur-gear  transmission  model  is  developed.  The  model 
includes  inertias  of  load  and  power  source,  stiffness  of  shaft,  time-varying 
mesh  stiffness,  and  damping  and  friction  inside  gear  transmissions. 

Governing  equations  of  the  model  are  derived  and  a  linearized  iterative 
procedure  for  solution  is  presented.  Parameter  study  including  rotating  1 

speed,  diametral  pitch,  applied  load,  damping,  stiffness,  and  inertia  will  be 
presented  in  NASA  TM- 100181  (AVSCOM  TM-87-C-3). 

•  , 
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